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Abstract 



5h . 

The effect of microwave irradiation on the phase dynamics of intrinsic Joseph- 
son junctions in high temperature superconductors is investigated. We demon- 
strate the influence of microwave's amplitude variation on the current- voltage 
characteristics and on the time dependence (temporal oscillations) of the 
electric charge in the superconducting layers. A remarkable changing of 
the longitudinal plasma wavelength at parametric resonance is shown. We 
demonstrate an effect of the microwave radiation on the width of the para- 
metric resonance region. 
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1. Introduction 

Intrinsic Josephson junctions (IJJs) in high temperature superconductors 
(HTSC) under microwave (MW) irradiation have attracted growing inter- 
est because of their rich physics and possible applications. IJJs are very 
complex, and their electrical and magnetic properties are strongly nonlinear. 
Hence a complete understanding of their physical dynamics is needed today. 
One of the most spectacular indications of the Josephson effect is locking 
of Josephson oscillations to the frequency of external MWs. As the Joseph- 
son frequency is proportional to voltage, the locking leads to appearance 
of steps at quantized voltages in the current- voltage characteristic (CVC), 
called Shapiro steps (SS) |2|. These steps appear in the CVC at discrete 
voltage values V n = nhu/2e, where u is the frequency of the applied MW 
signal and n is integer number. 
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It was shown by Koyama and Tachiki [3] that the system of equations 
for capacitively coupled Josephson junctions has a solution corresponding to 
the longitudinal plasma wave (LPW) propagating along the c-axis. So, the 
Josephson oscillations can excite the LPW by their periodical actions. The 
frequency of Josephson oscillations uj is determined by the voltage value 
in the junction, and at uj = 2ulpwi where ulpw is LPW frequency, the 
parametric resonance (PR) is realized. 

In this article we investigate the influence of external radiation on the 
parametric resonance of IJJs by: i) decrease of frequency of radiation which 
brings SS near PR region (PRR); ii) increase of amplitude of radiation. 
Bringing SS near PRR causes the width of this region to be shrink. An 
increase in the amplitude of radiation leads to (a) changing the wavelength 
of LPWs along the stack, (b) appearance of additional PR above SS. So, the 
external radiation essentially changes the physical picture of the coupled JJs. 



2. Model and Methods 

To investigate the phase dynamics of IJJ we use the one-dimensional 
CCJJ+DC model with the gauge- invariant phase differences <fi(t) between 
S-layers I and / + 1 in the presence of electromagnetic irradiation described 
by the system of equations: 

Vt - a(V l+1 + Vt-i - 2Vi) 

d l (1) 
/ - sin ipi - f3—±- + A sin ut + I noise 
at 

where t is dimensionless time, normalized to the inverse plasma frequency 
Up 1 , u p = a/ 2el c / HC, f3 = 1/a//3^, (3 c is McCumber parameter, a gives the 
coupling between junctions [3j, A is the amplitude of the radiation. To find 
the IVC of the stack of IJJ we solve the system of nonlinear second-order 
differential equations (1) using the fourth order Runge-Kutta method. In our 
simulations we measure the voltage in units of Vb = fou p /(2e), the frequency 
in units of u p , the bias current / and the amplitude of radiation A in units 
of I c . 

The time dependence of the electric charge in the S-layers is studied 
by Maxwell equation div(ee E) = Q, where e and e are relative dielectric 
and electric constants, respectively. The charge density Qi in the S-layer 
/ is proportional to the difference between the voltages V\ and Vi + \ in the 
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neighbor Josephson junctions Q\ = Q a(Vi + i — V{), where Q = es Vo/r^). 
To see clearly the effect of radiation only, we have taken even number of 
junctions in the stack N = 10 to escape the additional modulations of the 
electric charge oscillations in the PRR which appear in case of odd number of 
Josephson junctions. Numerical calculations have been done for a stack with 
the coupling parameter a — 1, dissipation parameter /3 = 0.2 and periodic 
boundary conditions. We note that the qualitative results are not sensitive to 
these parameter values and boundary conditions. In our simulations we put 
T m = 1000 for time domain, 5t = 0.05 for step in time, 51 = 0.0001 for step 
in current and we add to the bias current a small noise with amplitude in the 
interval (+10 -8 , — 10~ 8 ). The details concerning the numerical procedure are 
given in Ref. |4|. 

3. Results and Discussions 

It is known that in case of single JJ an increase in radiation amplitude 
A decreases a hysteresis region, i.e., it leads to the decrease of the critical 
current value and the increase of the return current Ir |5j. For the stack 
of coupled JJ the external radiation leads additionally to the series of novel 
effects related to the parametric resonance and the LPW propagating along 
the c axis 0, 0| . We will demonstrate below the changing of LPW wavelength 
and additional resonances around SS. In Fig. [T]we show the CVC for a stack 
with 10 JJs at a = 1 and /3 = 0.2 under radiation with amplitude A = 0.5 
and different frequencies. Curve 1 show CVC without irradiation. Inset (a) 
stresses the coincidence of the curves before SS. The SS does not appear at 
frequency smaller than u = 2.86, because before it a jump to another branch 
is happened. Inset (b) illustrates that, with decreasing the frequency of 
external radiation and approaching the parametric resonance condition, both 
the breakpoint current (hp) and the growing region (parametric resonance 
region) width are decreasing. Let us see the influence of the external radiation 
on the PR by the increase of amplitude of the radiation at fixed frequency. In 
this present work, we describe the case uj > 2ulpw, where the SS is above PR 
region in CVC. Fig. [2] shows 4 CVCs of a stack with 10 coupled JJs: without 
irradiation (curve 1) and under radiation with uj = 3, A = 0.3 (curve 2), 
A = 1 (curve 3) and A = 1.5 (curve 4). At u = 0, PR is characterized 
by breakpoint current hp — 0.5772, and breakpoint voltage Vb p — 28.65, 
corresponding to the Josephson frequency uj = 2.865. 

As shown in the figure, the first SS is developed on the outermost branch 



3 




Figure 1: CVC for the stack with 10 JJs under radiation with amplitude A = 0.5 and 
different frequencies. Curve 1 show CVC without irradiation. Inset (a) stresses the co- 
incidence of the curves before Shapiro step. Inset (b) demonstrates the changing of both 
breakpoint current (hp) and the growing region (GR) width with frequency variation. 



of CVC in the hysteresis region at V = ujj * N = 30. It was shown in Ref. 
that an increase of the amplitude of irradiation leads to the appearance 
of an additional PR before SS called radiation related parametric resonance 
(rrPR). Filled arrows indicate the positions corresponding to the appearance 
of a fundamental PR (fPR) in the stack, which realizes without radiation. 
Hollow arrow indicates the rrPR caused by irradiation. Inset (i) illustrate 
the modulation of the charge oscillations at A = 0.3, while inset (ii) shows 
the charge distribution among the S-layers in the growing region presented 
in inset (i). The charges on the neighbor layers are equal in magnitude and 
opposite in sign. This is true for all adjacent layers and corresponds to the 
7r-mode. So, the 7r-mode is survived under radiation with A = 0.3. Results of 
FFT analysis of the time dependence of voltage V(t) in each JJ and charge 
Q(t) in each S-layer show that this modulation is due to beating between 
the external and Josephson frequencies. We see in the inset (iii) an enlarged 
part of CVC corresponding to the appearance of fPR region in curve 3. 

The irradiation changes the shape of CVC as it is seen in the curve 3 
at uj = 3 and A = 1. At these parameters, as was mentioned above, the 
rrPR appears in the stack before SS (shown by hollow arrow). For this case, 
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Figure 2: CVC of the stack with 10 coupled J J without irradiation (curve 1) and under 
radiation with frequency u> = 3 and amplitude A — 0.3 (curve 2); A = 1 (curve 3); A = 1.5 
(curve 4). Inset (i) illustrate the modulation of the charge oscillation and inset (ii) shows 
the charge distribution among the layers at A = 0.3. Inset (iii) enlarges the fPR region in 
the curve 3. 



the charge-time dependencies in the fPR and rrPR regions are shown in Fig. 
3. Fig. 3(a) demonstrates the time dependence of the charge in the fPR 




Figure 3: Charge time-dependence at w = 3 and A = 1 for (a) fPR region; (b) rrPR 
region. The insets show the charge distributions among the S-layers in the growing region 
of resonances. 

region for the first S-layer at u = 3 and A = 1 which is combined with CVC 
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of the outermost branch. The filled squares mark the bias current steps in 
CVC. The inset shows the charge distribution among the layers, illustrates 
that the wavelength of the LPW is changed with increasing the amplitude 
to A = lOd. In Fig. 3(b) the time dependence of the charge in the rrPR 
region is introduced. We see clearly the correlation between the charge-time 
dependence and CVC. The inset shows that the wavelength of LPW at this 
rrPR is also A = lOd. 

4. Conclusions 

We demonstrated different effects of external electromagnetic radiation 
on the parametric resonance of coupled Josephson junctions. Approaching 
of the fundamental parametric resonance by the decreasing of the radiation's 
frequency shrinks the width of this resonance. An increase in the amplitude of 
radiation leads to (a) changing the wavelength of longitudinal plasma wave 
along the stack, (b) appearance of additional parametric resonance above 
Shapiro step. 
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